Crossover from Weak- to Strong- Coupling 
Superconductivity and to Normal State with 

Pseudogap. 



E. BabaevQ' 2 and H. Kleinert 1 

(1) Institut fur Theoretische Physik 
q\ , Freie Universitat Berlin, Arnimallee 14, 1000 Berlin 33, Germany 

ON \ (2) A.F.Ioffe Physico- Technical Institute 

Russian Academy of Sciences, Politechnicheskaja str. 26, St. Petersburg, 194021, Russia 

3 '■ 

Abstract 

\o ■ 

H ' For an electron gas with a <5-function attraction we investigate the crossover from 

weak-coupling to strong-coupling superconductivity as well as normal state near the 
\ temperature T*, at which the strong coupling produces a pseudogap in the energy 

spectrum due to the binding of electron pairs. We present curves for the behavior 
of the superconductive transition temperature, the gap formation temperature, the 

■^j- , gap size and several thermodynamic quantities as functions of coupling strength 

and temperature, both in two and three dimensions. 
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1 Introduction 

a; 

The crossover from BCS superconductors to a Bose-Einstein condensate of tightly bound 
fermion pairs was first studied many years ago in Refs. 0-0. Sparked by experimental 
j> . studies of short coherence length cuprate superconductors, it has recently attracted re- 
newed interest |4]]- []33|], the experiments show an anomalous behavior in the normal phase 
well above the superconductive transition T c (see the review [p4| , also |32|], [p5|1 -p8|). 



Anomalous is the temperature dependence of resistivity, specific heat, spin susceptibility, 
etc. . Moreover, angular-resolved photoemission spectroscopy experiments (ARPES ) in- 
dicate the existence of a pseudogap in the single-particle excitation spectrum |36| - [p9[] 
manifesting itself in a significant suppression of low-frequency spectral weight above T c , 
this being similar to the complete suppression below T c in an ordinary superconductor 
due to the energy gap. 



1 Emails: babaev@physik.fu-bcrlin.de, babaev@k-box.ioffe.rssi.ru; kleinert@physik.fu-berlin.de 



URL: http://www.physik.fu-berlin.de/~kleinertl Phone/Fax: 0049 30 8383034. 
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In essence, such anomalous properties of the normal state of superconductor can be 
described by a simple model of superconductivity in which electrons are bound to pairs 
by a 5-function potential. We merely must leave the BCS regime and go to stronger 
couplings 0. 

In this paper we present a detailed study of the crossover in such a model. Physically, 
the most important distinctions between conventional (BCS) and strong-coupling (Bose- 
Einstein) regime lies in the fact that in the former only a small fraction of the conduction 
electrons is paired with the superffuid density involving all pairs, whereas in the latter 
practically all carriers are paired below a certain temperature T*, although not condensed, 
effect that results in deviation from the Fermi liquid behavior in the region between T* and 
T c . The temperature has to be lowered further below some critical temperature T c < T* to 
make these pairs condense and establish phase coherence, which leads to superconductive 
behavior. We shall neglect the coupling to the magnetic vector potential throughout the 
upcoming discussion, so that the phase coherence below T c can be of long range, unspoiled 
by the Meissner effect which would reduce the range to a finite penetration depth. 

The pseudogap behavior between T c and T* is characterized by short-range pair corre- 
lation functions. The common physical origin of the superconductive gap below T c and the 
pseudogap above T c observed in cuprates is suggested by the above-quoted ARPES data, 
which show that the two gaps have the same magnitude and wave vector dependence. 
Important experiments on the gap properties are: 



1. In experiments on YBCO [50], [5T|, a significant suppression of in-plane conductivity 
v a b{u) was observed at frequencies below 500 cm -1 beginning at temperatures much 
above T c . Experiments fOfl , [|!| on underdoped samples revealed deviations from the 
linear resistivity law. In particular, a a b(w = 0; T) increases slightly with decreasing 
T below a certain temperature. 



2. Specific heat experiments [|4| also clearly display pseudogap behavior much above 



2 It should be noted that phase diagram in the case of such complicated materials as High-T c cuprate 
superconductors contains large region where anomalous properties of normal state are governed by anti- 
ferromagnetic correlations. 
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3. NMR and neutrons observations in jEJ and ||6| show that below temperatures T 
much higher than T c , spin susceptibility starts decreasing. 



4. Experiments on optical conductivity [plJl , |52|1 and tunneling exhibit the opening 
of a pseudogap. A review of actual experimental data confirming the pseudogap 
behavior of the underdoped and optimally doped cuprates is given in f5l| . 



In the model to be investigated in this paper, the crossover from BCS-type to Bose- 
type superconductivity and pseudogap state will take place either by varying the coupling 
strength, or by decreasing the carrier density 

Analytic calculations will be performed using a crossover parameter xq- This pa- 
rameter is directly related to the ratio of chemical potential and gap function at zero 
temperature. It is a monotonous function of coupling strength and carrier density [to be 
seen in Figs. [I] and|2| in two and three dimensions, respectively]. It is also a direct measure 
for the scattering length a s of the electron-electron interaction, at zero temperature in 
three dimensions [to be seen in Eq. (0)], and to the binding energy of the electron-electron 
pairs in two dimensions [to be seen in Eq. ([201)1. 

Our analysis of the model start with a mean-field approximation to the collective pair 
field theory f53| . It is well known, that mean-field results are reliable at all temperatures 
for weak coupling strength, i.e. in the BCS regime. As this regime is approached from 
the strong-coupling side, the temperature T* where pairs are formed and the temperature 
T c where phase coherence sets in merge to the single BCS phase transition temperature 
T c . This merging will be described analytically in this paper in two as well as three 
dimensions. 

Mean-field results are also reliable at strong couplings if the temperatures are suffi- 
ciently small to suppress fluctuations [see the discussion in 

Apart from that, mean-field results at stronger couplings seem to indicate correctly the 
position of the temperature T* where pseudogap forms due to precursor pair formation ||. 
The model shows the binding of noncodencerd pairs by the appearance a nonzero complex 



gap function A(T). The precise temperature behavior near T* is certainly predicted 
wrongly as being oc (T* — T) 1 / 2 . This would suggests a second-order phase transition, 
whereas the experimental data show a smooth crossover phenomenon. Here the effect of 
fluctuations is too important to be calculable analytically. 

In this paper we shall focus attention upon the onset of long-range order when lowering 
the temperature down from T*. Due to the strength of the coupling, this regime lies 
outside the mean-field approximation. It appears, however, that by extracting the lowest 
gradient terms governing the Gaussian fluctuations around the mean-field we shall obtain 
sufficient information to study the onset of long-range order and its destruction at T c < T*. 

In three dimensions, this was first investigated in Ref. by summing particle-particle 
ladder diagrams which correspond to Gaussian fluctuations around the mean field, in the 
functional integral formalism it was studied in Ref. |J. 

As it was noted above, at strong couplings, the formation of a gap does not imply 
the existence of superconductive currents. Although the modulus of the gap field A(x) = 
| A(x)| is nonzero, the phase of A will initially fluctuate so violently that long-range order 
cannot be established. This phenomenon is well-known from XY-models [0(2) classical 
Heisenberg model], which describe the fluctuations of a field of two-dimensional unit 
vectors n = (cos 0, sin 0) . The fluctuations of the phase angle <fi prevent the existence of 
a long-range order above a certain temperature. 

Much of the previous work on the model with an attractive (^-potential at strong 
couplings was devoted to zero-temperature or to a study of the critical region near T c where 
phase coherence sets in. Apparently, there is no analytic work displaying a complete set of 
global properties of the pseudogap state above T c . There only exist some numeric results 
in Ref. p|J on the paramagnetic susceptibility in the pseudogap phase of two-dimensional 
superconductors, derived within the same model after adding a paramagnetic term to 
Hamiltonian 0. In the present paper, we shall derive as well analytically the behavior for all 
coupling strengths of the pseudogap as well as the thermodynamics of the pseudogapped 
phase in both two and three dimensions. 



For Monte-Carlo study of the normal state just above T c see Ref. 30 
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The weak- to strong-coupling crossover of the Kosterlitz-Thouless in two dimensions 
was studied in Refs. ||26|| , [23|. In Ref. it was investigated within the same model as 



ours at a fixed carrier density, but only numerically^ . The different properties of size and 
phase fluctuations was also exploited in Refs. ]54] and f^ . 
The plan of the paper is the following: 

We begin by reproducing the results of Refs. [O], ]T^| for the gap A(0) and chemical 
potential at zero temperature. These are subsequently extended by equations for the 
temperature behavior of gap and pseudogap, as well as of thermodynamic functions. 

In Section 2 we calculate the crossover from BCS to strong-coupling superconductivity 
at small but finite temperatures and present solutions for the gap and thermodynamic 
functions of the superconductive state. 

In Section 3 we study the crossover from BCS superconductivity near T c to the onset 
of pseudogap behavior near T*. 

In Sections 4 and 5 we go beyond mean-field approximation and study crossover of 
superconductive transition in both two and three dimensions. In the strong-coupling limit 
of superconductors, we set up an equivalent XY-mode\ in two as well as three dimensions 
with the help of a gradient expansion for the phase of the order parameter. This allows 
us to find the onset of phase coherence in a non-perturbative way. Even though in two 
dimensions, this was done before via Kosterlitz-Thouless arguments we argue that merging 
of temperatures of pair formation and condensation was missed. In three dimensions, our 
discussion is new. Previous results were derived form a study of the condensation of the 
gas of composite bosons via retaining corrections to the number equation. 

In the weak-coupling limit of both two- and three-dimensional supercondictors we 
show how the XF-model transition converges with the transition in the BCS theory. 



4 We shall see in Section 4 that these numerical results do not cover the entire crossover region, in 
particular, the merging of Tkt and T* in the weak-coupling region is missing - the effect that we show 
analyticaly in our paper. 
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2 Crossover from BCS to Strong-Coupling 
Superconductivity Near Zero Temperature 



The Hamiltonian of our model is the typical BCS Hamiltonian in D dimensions (H — 1) 
H = J2 Jd D x^(x) [~ -(A M*)+9 /rf D x^(x)^(x)^(x)^ T (x), (1) 

where Vv( x ) is the Fermi field operator, a =|, \ denotes the spin components, m is the 
effective fermionic mass, and g < the strength of an attractive potential g5(x — x'). 

The mean-field equations for the gap parameter A and the chemical potential \i are 
obtained in the standard way from the equations (see for example [53], and Ap- 
pendix A): 



- = -V — tanh^, (2) 



n - I ? (l-|t«h|], (3) 



where the sum runs over all wave vectors k, N is the total number of fermions, V the 
volume of the system, and 



£ k = ^k + A 2 with (4) 

are the energies of single-particle excitations. 

Changing the sum over k to an integral over £ and over the directions of k, on which 
the integrand does not depend, we arrive in three dimensions at the gap equation: 

- = K 3 f d£ — ^ = tanh — ' ~j~ - , (5) 
g J-» 2V£ 2 + A 2 2T ' K ' 

where the constant n% = m 3 / 2 /^/2ir 2 has dimension energy -3 / 2 /volume. In two- 
dimensions, the density of states is constant, and the gap equation becomes 



- — 1^2 ( d£ — = tanh — - 1" , (6) 

g J-im 2V£ 2 + A 2 2T 
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with a constant n 2 = m/27r of dimension energy 1 /two- volume. In two dimensions, the 
particle number in Eq. (0) can be integrated with the result: 



m 



n = — {J fx 2 + A 2 + fi + 2Tlog 



1 + exp 



+ A 2 



(7) 



the right-hand side being a function n(fx, T, A). 

The 5-function potential produces an artificial divergence and requires regularization. 
A BCS superconductor possesses a natural cutoff supplied by the Debye frequency uj d . 
For the crossover problem to be treated here this is no longer a useful quantity, since 
in the strong-coupling limit all fermions participate in the interaction, not only those in 
a thin shell of width ud around the Fermi surface. To be applicable in this regime, we 
renormalize the gap equation in three dimensions with the help of the experimentally 
observable s-wave scattering length a s , for which the low-energy limit of the two-body 
scattering process gives an equally divergent expression P-|T3[: 



m 1 1 ^ m 

4^ ~ 3 + 



Eliminating g from (J8|) and (Q) we obtain a renormalized gap equation 

(9) 



"' 1 £ 



Ana s V k 



1 E k m 
tanh 



2E k 2T k 2 . 

in which 1 / kpa s plays the role of a dimensionless coupling constant which monotonically 
increases from — oo to oo as the bare coupling constant g runs from small (BCS limit) 
to large values (BE limit). This equation is to be solved simultaneously with (Q). These 
mean-field equations were first analyzed at a fixed carrier density in Refs. and [0] . Here 
we shall first reproduce the obtained estimates for T* and fi. 

In the BCS limit, the chemical potential fi does not differ much from the Fermi energy 
e_F, whereas with increasing interaction strength, the distribution function broadens 
and /j, decreases, and in the BE limit, on the other hand we have tightly bound pairs 
and nondegenerate fermions with a large negative chemical potential 3> T. Analyzing 
Eqns. (HD and @ we have from (|3|) for the critical temperature in the BCS limit (fx ^> T c ) 
t bcs = 8e- 2 e 7 7r- 1 e F exp(-7r/2A; F |a s |) where 7 = -r'(l)/r(l) = 0.577. . . , from (|) we 



have that chemical potential in this case is /i = e F . In the strong Eq. (pj) determines T*, 
whereas Eq. (|3|) determines \i. From Eq . (^) we obtain that in the BE limit fi = — -Eft/2, 
where = l/ma 2 s is the binding energy of the bound pairs. In the BE limit, the 
pseudogap sets in at T* ~ Eb/2\og(Eb/e F ) 3 ^ 2 . A simple "chemical" equilibrium estimate 
(fib = 2/i/) yields for the temperature of pair dissociation: Tdissoc — Ebj \og(Eb/e F ) 3 ^ 2 
which shows at strong couplings T* is indeed related to pair formation M, || (which in 
the strong-coupling regime lies above the temperature of phase coherence pfl-|32]|). 

The gap in the spectrum of single-particle excitations has a special feature ||, |7| 
when the chemical potential changes its sign. The sign change occurs at the minimum of 
the Bogoliubov quasiparticle energy E^ where this energy defines the gap energy in the 
quasiparticle spectrum: 

^2 , A 2\V2 



E gap = min (£k + A J . (10) 

Thus, for positive chemical potential, the gap energy is given directly by the gap function 
A, whereas for negative chemical potential, it is larger than that: 

_J A for /i>0, . . 

^P- 1(^2 +A 2) 1/2 for /i<Q _ W 

In three dimensions at T = 0, equations (pj), (H) were solved analytically in entire crossover 



region in [11 to obtain A and fi as functions of crossover parameter l/k F a s . The results 
are 



A 



e F [x h(x ) + I 2 (x )] 2 / r ' 

jJL fl A Xq 



e F A e F [xohixo) + / 2 (xo)] 2/3 ' 
1 4 x q I 2 (xq) - h{x Q ) 



k F a s 7T [x Ii(x ) + I 2 (x )] 1/3 '' 



(12) 
(13) 
(14) 



with the functions 



POO 

h(xo) = dx 
Jo 



x 2 



(X 4 - 2X X 2 + X 2 Q + l) 3 / 2 



1 + xD^eC-, «) - - 2n ) 2 , 1/4 F(-, k) 

V 2 ' Ax\{l+xiyi A 2 y 




Figure 1: Gap function A and chemical potential /i at zero temperature as functions of 
xq in three dimensions. 

In two dimensions, a nonzero bound state energy eo exists for any coupling strength. 
The cutoff can therefore be eliminated by subtracting from the two-dimensional zero- 
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temperature gap equation 



1 1 1 m f°° , 1 , . 



where 2; = k 2 /2mA — x , the bound-state equation 



00 



1 1 ~v 1 m I d x (20) 



After performing the elementary integrals, we find: 



| = v/l + xg-xo. (21) 

From Eq. ([?]) we see that at zero temperature, gap and chemical potential are related to 
%o by 

A 2 

r = — 7r=i' (22) 

^ = (23) 

€f Xo + V 1 + ^0 

The two relations are plotted in Fig. ^. Combining (|2~lD with (^) we find the dependence 
of the ratio eo/e_F on the crossover parameter xq\ 

eo 2 V^+g-^o (24) 



€F V 1 + x o + x ° 

We have extended all these relations to non-zero temperature. For this purpose, we 
do not fix the carrier density but assume the presence of a reservoir which provides us 
with a temperature-independent chemical potential \i = fi(l/kpa s ;T = 0). Such a fixed 
\x will be most convenient for deriving simple analytic results for the finite-temperature 
behavior of the system.^ 

5 In Ref. J24|] , the temperature dependence of the chemical potential was calculated numerically within 
a "fixed carrier density model", where it turned out to be very small in comparison with the dependence 
on the coupling strength. 
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Figure 2: Gap function A and chemical potential /i at zero temperature as functions of 
Xq in two dimensions. 



In our calculation we use xq as the most convenient crossover parameter, since it 
depends via the simple relation (fl8|) on the chemical potential which can be measured 
rather directly experimentally [EJJ. The parameter xq ranges from — oo in the strong- 
coupling (Bose-Einstein) limit to oo in the weak-coupling (BCS) limit. The relation 
between x and the inverse reduced coupling strength between the electrons l/kpa s is 
plotted for three-dimensional system in Fig. ||]. The corresponding relation ( p4|) in two 
dimensions between x and the bound state energy e of the electron pairs is plotted on 
Fig. |. 

In Fig. |5] shows the temperature behavior of A near T = for different coupling 
strengths in three dimensions. In Fig. [£] does the same thing in two dimensions. Figures [7] 
and H display dependence of the temperature T* where the gap vanishes on the coupling 
strength parameter x . 



11 




-60 -40 -20 20 T 40 60 



Figure 3: Dependence of l/kpa s on the crossover parameter x in two dimensions. 



As was noted above, recent experiments on the underdoped cuprates showed that an 
ordinary superconductive gap develops smoothly to a pseudogap above T c . Within the 
mean-field approximation we obtain the analytic results shown in Fig. |5| and ^| 

Let us now turn to the region near zero temperature, where we can derive exact results 
for gap. From flT2] ) we extract the asymptotic behavior in the three-dimensional case for 
xq > 1. In this region one can assume density of states to be roughly constant, since the 
integrand of (|5]) is peaked in the narrow region near £ = 0. The small-T behavior is 



A(T) = A(0)-A(0) A /- 





~ A(0)1 


/A(0) eXP 


T 



1 + erf 



V 



arg + 1 - 1 
\| T/A(0) 



\ 



(25) 



where erf (x) is the error function. Since the density of states is nearly constant in this 
limit, the same equation holds in two- dimensions — apart from a modified gap A(0) given 
by©. 

In the weak-coupling limit, Xq = /i/A(0) tends to infinity, and the expression above 
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Figure 4: Dependence of eo/e^ on the crossover parameter xq 



approaches exponentially fast the well-known BCS result: 



A(T) = A(0) - [27rA(0)T] 1/2 exp[ 



A(0) 1 

T 1 



(26) 



For strong couplings with Xq < — 1, the three-dimensional integrands are no longer peaked 
in the narrow region so that the density of states can no longer be taken to be constant. 
Taking this into account, we find: 



A(T) = A(0) - ±y/=^ (^) 3/2 exp 



fi 2 + A 2 (0) 



T 



(27) 



From Eq. ([27]) we see that near T = the gap A(T) tends in the strong-coupling limit 
exponentially to A(0), forming plateau near T = 0. 

In two dimensions we arrive at similar result: an exponentially growing plateau near 
T = in the strong coupling limit: 



A(0) / \ A(0) 2 + /i 2 
A(T) = A(0) - -U^ ( V K, T 

where E\ is the exponential integral E\(z) = J z °° e^/t dt. 



{21 
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0.4 , 0.6 

Figure 5: Temperature dependence of gap function in three dimensions. Solid line 
corresponds to crossover parameter x = 10 (i.e., in the BCS regime), the crosses to 
xq = (i.e., in the intermediate regime), lines with boxes and circles represent xq = —2 
and Xq = —5 cases correspondingly and the dashed line corresponds to x = — 10 (i.e., in 
strong-coupling regime). 



For very strong couplings, Eq. becomes: 

A(T) = A(0) - ^> , T exp 



/i 2 + A(0) 2 



(29) 



2 ^/p 2 + A 2 (0) V T 

Let us also calculate thermodynamical quantities near T = 0. For the thermodynamic 
Gibbs potential fi(T, /i, V) we calculate 



Here and in the sequel in this section, A(0) will be replaced by A. In three dimensions, 
Eq. ( |30"1) turns into the 

A 2 



A 2 J& + A 2 

^=^= tanh 2T log 

2T & 



2 cosh 



2T 



(30) 



n 
v 



K 3 I d^J^ + fi 

J—u 



2ve + a 2 



+ A 2 / v^ 2 + A 2 \ 
tanh ^± 2Tlog 2 cosh + £ 



2T 



2T 



(31) 
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Figure 6: Temperature dependence of gap function in two dimensions. Solid line corre- 
sponds to crossover parameter x = 10 (i.e., in the BCS regime), the crosses to xq = 
(i.e., in the intermediate regime), lines with boxes and circles represent xq = —2 and 
x = —5 cases correspondingly and the dashed line corresponds to x = —10 (i.e., in 
strong-coupling regime). 



In two dimensions, we obtain instead: 



v = K2 



A 2 JP + A 2 / Jf 2 + A 2 \ 

^tanh^± 2Tlog 2cosh^± ] + £ 



(32) 



2y/C 2 + A 2 IT "° V IT 

We regularize the thermodynamic potential Q s of the condensate subtracting Q n = Q(A = 
0). At T = and for weak couplings this is found to depend on temperature as follows: 



V 



v 



n , — 
- = K3VA* 



A2 1 , 1 n — iTo 



In the BCS limit (x — > oo) this reduces to the well-known result 

n s r a 2 " 

In two dimensions, we have a formula valid for any strength of coupling: 



v 



= K 2 



\^^ 2 + A 2 



(33) 



(34) 



(35) 
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Figure 7: Dependence of T* on crossover parameter in three dimensions. 



with the BCS limit 



V 



K-2 



A2 

2 



(36) 



In both three- and two-dimensional cases in the BCS limit we can write a small temper- 
ature correction ttT 2 /3 to the thermodynamic potential. In the opposite limit of strong 
couplings, we find in three dimensions the strong-coupling limit: 



v = -^ A5/2 (-*») 



-3/2 



(37) 



The gap A(0) has by Eq. (11) the strong-coupling limit A(0) ~ ei?[16/37r] 3 / 2 |xo| 1 ^ 3 yielding 
the large-Xo behavior 



V^_ k 5/2 *_( 16 \ 15/4 

V ~ Ks€p 64 Utt; 



-2/3 



l-n,| - (38) 

In two dimensions, we substitute the gap function A of Eq. (^) , into the thermodynamic 
potential Q55D, and obtain for strong couplings where fj, < 0: 

n 



v 



= o 



(39) 
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Figure 8: Dependence of T* on crossover parameter in two dimensions. 

Let us now turn to the entropy. In three dimensions near T = it is given for weak 
couplings by: 

Y 



s_ 

V 



, 2ttA3 / A 





( 


1 + erf 









\ T/A 



J 



+ 2/^exp 



T 



,(40) 



For /i/A — > oo, this reduces correctly to the BCS result: 

S ^ /8ttA 3 / A 



exp ( - — 



(41) 



In two dimensions, the result is similar [with A = A(0) given by Eq. 



s 

— = 

V 



'27rA 2 



A 







1 + erf 









^0 + 



\ T/A 



J 



2/iexp 



T 



,(42) 



The BCS limit of this is 



v = K2] 



/8ttA 3 / A 
— exp (-- 



(43) 



17 



In the strong-coupling limit where ///A <C — 1, we have for the entropy in three dimensions 

(44) 



S I 2 , A2 I vV + A 2 * 
y = K 3—T 1 ^H 2 + A 2 exp I 



and in two dimensions: 



S ( y//F+A^ 
— = -2/t 2 /i exp I 



(45) 



From the entropy, we easily derive the heat capacity at a constant volume cy. In three 
dimensions it is given near T = for weak couplings by 



Cy = K 3v //iV / 27rA3 < 



7 -' 26XP { t ) 



1 + erf 



xl + 1 - 1 



^ T/A 



/J 



(46) 



reducing in the limit Xo — > oo to the BCS result 

cy = k 3 ^J2V2ttA 3 ^ exp (-^) . 
In two dimensions, the weak-coupling behavior is 

A 



(47) 



Cy = K 2 V2nA 3 < 



7 2 rXP { T 



1 + erf 



xt + 1 - 1 



\| T/A 

while the strong-coupling behavior in three dimensions is 

c,^3^T-V 2 (, 2 + A 2 )exp(-^±^) 



l\) 



(48) 



(49) 



and in two dimensions 



cy = 2k 2 — exp I 



(50) 



3 Crossover from BCS Superconductivity Near T c to 
Onset of Pseudogap Behavior 

We now turn to the region near T*, for which we derive asymptotic behavior of the ratios 
A(T)/T* and A(T)/A(0) as well as thermodynamic quantities. In doing so, we shall 
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consider A(T)/T as a small parameter of the problem. In the calculations near T* it is 
convenient to use fi/2T* (this ratio tend to oo in the weak-coupling limit and to the — oo 
in the strong coupling one) as a crossover parameter rather than Xq. In three dimensions, 
we find for weak couplings 



ACQ 
2T* 



T 



1 + tanh 



2T* 



tanh 



In the limit (jt/2T* 



ji/2T* (fi/2T*) 2 2T 
-> oo this tends to the BCS result 
A(T) / T 



2 \ 2 
7T, 



, 2 U 

1 H — arctan 

7T 7rT" 



(51) 



(52) 



In the opposite limit of strong couplings, T* and A(0) tend to infinity. The ratio A(T)/T 
near T* tends to zero exponentially as a function of the crossover parameter /i/2T*: 

16 r. t\, u \ 3 / 2 (53) 



A(T) 



2T* 



'2ir 



T 



2T* 



In two-dimensions the near-T* formula (pTD holds over the crossover region: In the weak- 
coupling limit, this formula reproduces the BCS result Q52p. In the strong-coupling limit, 
we find as in three dimensions a ratio A(T)/T which tends to zero exponentially as a 
function of the crossover parameter fi/T* : 



A(T) 



2T* 



7T/ 



V2T* 



1 - £ ) e"' T 
T* 



Let us calculate the dependence of T* on the crossover parameter /i / 2T* in the strong- 
coupling limit. In three dimensions, we obtain from Eq. (|5) the relation 



]_v 2/3 



exp 



This is solved for T* (up to a logarithm) by 



^---/ilog" 1 



2 /! 

3T* 



.if. 



(55) 



(56) 
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[see also the discussion after the formula @]. As a function of the crossover parameter 
xq we obtain 

T* ^ 1 /16\ 2 / 3 
7^~2 V3W 

In two dimensions we find from (0) 



(^^Vol^log-^V^kol). (57) 



e F 2 \ T 

and thus 



7 'expf-^l. (tt 



T^-zilog-^-A) (59) 
As a function of Xo, this implies 

— = 2x 2 Q !og _1 (2V2|x |) . (60) 

Let us also derive in the strong-coupling region the dependence of the ratio A(0)/T* on 
the crossover parameter which in three dimensions reads 

A<0) - 4 < ^W^, (61) 



T* v^r V TV \2T* J ' 

and in two dimensions: 



A(0) -4f-^V /! exp^). (62) 



In the weak-coupling regime, both three- and two- dimensional cases yield the result 

The temperature T* is in the weak-coupling regime of three- and two dimensional systems 
the following function of x$: 

Tl^l(L-±\ (64 ) 

Using these results, we can also calculate the asymptotic behavior of the ratio 
A(T)/A(0) near T*. In three dimensions, the strong-coupling limit yields 



A(T) 
A(0) 



2 



/ • (65) 



2 V 2T* / V T* 
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In two dimensions we have that in the strong-coupling limit this ratio tends to 

2 , , „ -, . _1 , rp 

1- 



A(T) 
A(0) 



1/4 1 
8 V^ 2 ~ 4 



(66) 



At weak couplings, both three- and two-dimensional gap functions are given by 



A(T) 
A(0) 



4vr 2 
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1 + tanh 
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tanh 



, 2 u 
1 H — arctan — — 

7T 7TT* 



(67) 



/i/2T* W2T*) 2 2T* 
In order to calculate thermodynamic potential near T* we expand the general expression 
by A(T)/A(0) and keeping the terms of the same order we get : 

J (2tt) d 2T* 

2 e 



V 



!/ 7 (2-)' 
d^p 1 / 1 



A* r d D p 1 / 



cosh 



- tanh — |- 
£ 2T* 



(68) 



2T* 2T* 

where -D is the space dimension. Note once more that we consider here the temperature 
evolution of the system at a fixed chemical potential (/i(T, l/k F a s ) = /x(0, l/&Fa s )) and 
regularize f2 by a subtraction of fi n = Q(A = 0). 

The thermodynamic potential in three dimensions in the weak-coupling regime near 
T* is given by 



ft. 
V 



(T* - T)A 2 
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A 4 
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2T* 
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1 + tanh 
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tanh 



(69) 



4 \/i/2T* (/i/2T*) 2 2T*y v. 
In the BCS limit, this reduces to the well-known formula: 

ft* 

V 

In the strong-coupling limit we have: 

H 



2 u 

— arctan 

tt irT* 



r-.l( T 1 A 2 ^ 



y = -a ^A 4 (2TT 3/2 



2T* 
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3/2 / rp s 



(70) 



(71) 



Using the asymptotic estimates derived above for the strong-coupling limit, and the fact 
that in this limit 



JL 

rp* 



log 



.JL 



-const x log(jx |), 



(72) 
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we find near T* the difference between the thermodynamic potential of the gapless and 
pseudogapped normal states: 



a 



. — —const 1 
V V T* 



T 



\x Q 



-3/2 



(73) 



In two dimensions near T = T*, the thermodynamic potential of the gas of pairs is given 
by the formula holding for the crossover region 



v = ~ Kr 



(T* - T)A 2 
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1 + tanh 
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(74) 



+ 



A 4 



4 (2T 
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tanh 



4 \fi/2T* (/u/2T*) 2 2T*/ v zr 
In the BCS limit, this yields the familiar result 

T 1 A 2 



2 u 
— arctan ■ 



7T 



TFT* 



-K2A 1 



and in the strong-coupling limit: 

a f r 



^ /, 2 j^l--jA 2 exp Vr ^ , 4 (2 ^ )2 

Using the earlier-derived asymptotic behavior plus the limiting equation (|72|) which also 
holds for two-dimensional case, we derive for the thermodynamic potential the xo-behavior 



T* 2vr 2 T 



A 4 



2 T* 2 / ' 



(75) 



4 7T 2 / /i/2T 



(76) 



a /. v > 



— ~ -const x 1 1 - -1 log |x |. 
The entropy behaves near T* in three dimensions in the weak-coupling regime like 



(77) 



S s S S n 
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r A 2 
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( V 



\2T* 
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with the BCS limit 



S s ^A 2 

v = ~ k ^y;: 



(79) 



The opposite strong-coupling limit is in three dimensions: 

\^2rp*-X/2 exp ^ 



Ss 

V 



-k 3 - 
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Inserting the above asymptotic formulas for A, //, T*, we find 



V 



t4 ~ -const x I 1 — — I \x \ 5/3 • 



(81) 



In two dimensions, the entropy is given in the entire crossover region by 



,2 r 



= — K*2 

V 2T* 



and has the BCS limit 



1 + tanh 



A 2 



JL 

IT* 



(82) 



V ~ K2 T*' 



(83) 



while the strong-coupling limit becomes: 



9 A 2 
V T* 



(84) 



Using corresponding asymptotic formulas for A, /i, T* in two dimensions, this depends 
on xq as 



tt = -consl x I 1 — — ] x 



V 



'0 



(85) 



In order to derive the specific heat we must take into account the temperature dependence 
of the gap. 

In three dimensions, we find in the weak-coupling region near T*: 

2 
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— = 2Tk 3 t/J1j 



1 + tanh 



JL 

2T* 



1 



tanh ■ ^ 



H/TT* (n/2T*) 2 2T 
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1-1 — arctan — — - 
n ttT* 



(86) 



which has the well-known BCS limit: 

^ ~ K 3v fji7T 2 T c . 

In the strong-coupling limit, we find in three dimensions 
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Inserting earlier derived asymptotic formulas we see that C s tends in the strong-coupling 
limit to zero like 

C 

-^r ~ const x |x |~ 2 (89) 
In two dimensions, the result for the entire crossover region reads 



(~i 1 + tanh 

" 2T*k 2 



2T* 



V 1 



1 1 , u 

tanh ■ 



H/2T* (fi/2T*) 2 2T* 
This becomes in the BCS limit 

C, 



. 2\ / 2 fi 
+ — 1 H — arctan ■ 



71 J \ 7T 7rT* 



(90) 



r ^ 2 vr 2 T*, (91) 



and in the strong-coupling limit 

C s (I 4X- 1 (2ft 



^(--- 2 ) exp(£). (92) 



As a function of x , the result is 

C 

— * ~ const x Xq 2 . (93) 

From the above calculation near T* we see that both quantities S s and C s tend quickly 
to zero with growing coupling strength in the pseudogapped regime (like a power of the 
crossover parameter \xq\ or exponentially as a function of crossover parameter fi/2T*). So, 
at very strong couplings T* is getting less and less pronounced with increasing coupling 
strength. 



Note that in the strong-coupling regime, the modified gap function fi 2 + A 2 [see 
Eq. |TIJ enters the expressions for thermodynamical quantities below T* the same way as 
an ordinary gap in BCS limit [see Eqs. (II), (HH), (M), (EDI)]. 



4 Phase fluctuations in Two Dimensions 
and Kosterlitz-Thouless Transition 



In the previous sections we have calculated the properties of the model in the mean-field 
approximation. 
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Now we are ready to go beyond the mean-field approximation, which as discussed 
in the introduction, supply us with modulus of the gap function. In this chapter we 
make use of derivative expansion which determines the crucial stiffness parameter for 
the study of phase fluctuations, that in two dimensions leads to the Kosterlitz-Thouless 
transition, at which the expectation of the complex order field A(x) = |A|(x)e l0 ( x ) van- 
ishes in the pseudogap state. In these calculations we assume with other authors (]2"4|. 
|54|| ) that the phase fluctuations do not significantly affect modulus of A. We shall first 



study the two-dimensional system, where the mean-field solution receives the strongest 
modifications from the violent phase fluctuations, as articulated by the Coleman- Mermin- 
Wagner-Hohenberg theorem 1)7]] which forbids the existence of a strict long-range order, 
leading to a power behavior of correlation functions for all temperatures below Tkt- 

The crossover of the Kosterlitz-Thouless transition from weak to strong coupling was 
first considered in [|TJ|, 2E|, and studied recently by means of an XY-model in with 
the stiffness derived from a fixed nonvanishing modulus of the order parameter A. In 



Appendix A, we outline the derivation of the effective Hamiltonian in [24]. Under the same 
assumptions, we shall analyze the weak- to strong-coupling crossover of the Kosterlitz- 
Thouless transition in the present work. 

Writing the spacetime-dependent order parameter as A(x)e l0 ^ x ', where x denotes the 
four- vector x = (r, x) formed from imaginary time and position vector, the partition 



function may be written as a functional integral |53], |6(| |U 

Z(p,T) = J AVAVOexp T,A(x),d9(x))\, 



where 



pn(n,T,A(x),d6(x)) = - [ dr [ dxA 2 (x) - Tr logtT 1 + Tr logGo 1 

g Jo J 



(94) 



(95) 



is the one-loop effective action, containing the inverse Green function of the fermions in 
the collective pair field 



T 3 



nA(r, x) 

(W(r,x)) 2 l f riV 2 fl(r,x) zW(r,x)V 



+ 1 



4m 



2m 



(96) 
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Here r 1; r 3 are the usual Pauli matrices, and G = G| Mi A,e=o is added for regularization. 

Let us now assume this smoothness implying that phase gradients are small. Then 
Q(fi,T, A(x),d9(x)) can be approximated as follows: 



Q(fM, A(x), dd{x)) ~ 0,^0*, T, A, 00(e)) + ClpotO", T, A), 



(97) 



A=const 



(9? 



(99) 



with the "kinetic" term (see @, @) 

OO 1 

and the "potential" term 

Q pot (fi, T, A) = Q | d D xA 2 - TTr log^T 1 + TTr logG ^ 

The latter coincides with our earlier mean-field energy (see also Appendix A), determining 
the modulus of A(/i, T) and thus the stiffness of phase fluctuations. The kinetic part Jl kin 
contains gradient terms whose size is determined by the modulus of A(/i, T). Given the 
stiffness, one may immediately set up an equivalent XY-model. Both £\m and f2 pot are 
expressed in terms of the Green function of the fermions, which solves the equation 



A=const 



- I9t + T3 _ + „ +T1 A 



and 



E{d6) = r 3 



~id T e (V0) 2 " 

2 ' 8m 
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0(r,x) = <J(r)5(x) 
iV 2 iV0(r,x)V 
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The gradient expansion that we use to determine stiffness was first made in Ref. |60| at 
zero temperature. In Ref. [^|, the kinetic term f2kin was calculated at finite temperature 



for arbitrary chemical potential retaining terms with n = 1, 2 in the expansion fl9~8 
The result is (see Appendix A) 
T rP 



fikin = -J q drj d D x [n(/i,T,A)z9 T + J( /U ,T,A(/i,T))(V0) 2 + K(/i,T,A(/i,T))(9 r 0) 2 

(102) 

where J(/i, T, A) is the stiffness coefficient 



J(/i, T, A) = ^-n(ji, T, A) - /°° x + »/ 2T 

4m VP ' 47rJ- M /2T cosh 2 Jx 2 + A 2 /4T 2 



(103) 
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^,r,A) = ^(l + 7 =^ta„h^!), (104) 

and n(n, T, A) is the density of fermions (0) which varies with temperature in our model. 
At the temperature T* where the modulus of A vanishes, also the stiffness disappears. 
The kinetic term corresponds to an XY-model with a Hamiltonian [BI|, |J: 

H= J -j dx[W(x)] 2 , (105) 

the only difference with the standard XY-model lying in the dependence of the stiffness 
constant J on the temperature, which is determined from the solutions of gap and number 
equations (§) and (|3|). Clearly, in this model the Kosterlitz-Thouless transition always 
take place below T*. In the XY-model with vortices of a high fugacity, the temperature 
of the phase transition is determined by a simple formula p2" |: 



P KT = \j (106) 



which follows from the divergence of the average square size of a vortex-antivortex pair. 
Since these attract each other by a Coulomb potential v(r) = 2n Jlog(r/r ), the average 
square distance is 

roc 1 

<r 2 >oc/ drrr 2 e- (2 ^ T ^ r ^ oc - — , (107) 

Jr 4 — 27TJ/T 

which diverges indeed at the temperature ( |106j ). In our case T KT should be determined 
self-consistently: 

T KT = \ J{fi, Tkt, A(/i, T kt )). (108) 

From ( |103| ), (|^) and ( |106| ) it is easily seen that Tkt indeed tends to zero when the pair 
attraction vanishes in which case A(T = 0) = 0. In general, the behavior of Tkt for 
strong and weak couplings is found by the following considerations. We observe from 
the above-derived limiting formulas for A(T, jj) and T*, that the particle number n does 
not vary appreciably in these limits with temperature in the range < T < T*, so that 
weak-coupling estimates for Tkt derived within the model with temperature-independent 
chemical potential (i.e. when the system is coupled to a large reservoir, see discussion after 
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the formula fl23|) ) practically coincide with those derived from a fixed fermion density. 
Further it is immediately realized from the equations ( |108| ), (|?p, ( |103| ) and (|J) that in the 
weak-coupling limit A(Tkt, /-O/Tkt is a small parameter. At zero coupling, the stiffness 
J(/i, Tkt, Tkt)) vanishes identically, such that an estimate of J at weak couplings 
requires calculating a lowest-order correction to the second term of eq. (|103|) proportional 
to A (Tkt, /^) /Tkt- Thus weak-coupling expression for stiffness reads [ 031 : 

7C(3) A(T) 2 



Using the limiting behavior (51) of A(T) [see also discussion after fl53"|) 1 we find after some 
algebra the weak-coupling equation for T K t: 



e F /„ T, 



KT 



where €f = (ir/m)n is the Fermi energy of free fermions. It is useful to introduce reduced 
dimensionless temperatures Tkt = Tkt/^f and T* = T* /ep which are small in the weak- 
coupling limit. Then we rewrite Eq. ( |110| ) as 



KT 



41 + 1/4T* 
For small T* we may expand 



(111) 



f KT ~ T* - (112) 

This equation shows nicely how for decreasing coupling strength Tkt merges with T*. 
As a function of the crossover parameter xq, the temperature Tkt behaves like 

e 7 1 

T KT « • (113) 

7T X 

The merging of the two temperatures in the weak-coupling regime is displayed in Fig. |^. 

Consider now the opposite limit of strong couplings. There Eqs. ( |108j ), (0), ( |103| ), 
and (H) for Tkt, n(T, /i), and A(T, /x) show that Tkt tends to a constant value. From 
Eqs. (TO, (|), (0), and (TO, as well as the limiting expressions ([28]) and (^) it follows 
that in the strong-coupling limit A(Tkt) is always situated close to the zero-temperature 
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Figure 9: Weak-coupling behavior of T KT (x ). The solid line is T*, the dashed line 
represent Tkt- 

value of A (T K t, A*) ~ A(T = 0, //). Taking this into the account we derive an estimate for 
the second term in ( |1U3| ), thus obtaining the strong-coupling equation for Tkt: 



7i I 1 T KT 

T K t ^ - < — n exp 

8 1m 7r 

which may further be expanded as 

7r J 1 T KT 
Tkt - - < — n exp 

8 m 7T 



/i 2 + A 2 (T KT ,/i) 



Ti 



KT 



/i + A 2 (T KT ,/i) 



(114) 



(115) 



Tkt 2/iT K T 

With the approximation A(T K t,/-0 ~ A(T = 0,/i) and the limiting behavior ( p2|) we 
find that the first term in the exponent tends in the strong-coupling limit to a constant, 
A 2 (T K t, Ai)/2/iT K T — ► —4, whereas the first term in the brackets tends to — oo, so that 
Eq. ( |115| ) has the limiting form 



Ti 



7r n 



KT 



< 1 exp 



As a function of xq crossover parameter, this reads 



T 



KT 



7r n 
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1 exp 
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(116) 



117) 
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Thus for increasing coupling strength, i.e., decreasing crossover parameter xq <C — 1, the 
phase-decoherence temperature Tkt tends very quickly towards a constant: 



5 m 

In this limit we know from Eq. (|^) that the difference in the carrier density at zero 
temperature, n(T = 0), becomes equal to n(T = Tkt)-, so that our limiting result coincides 
with that obtained in the "fixed carrier density model": 

e F (no) 7T 

T KT = — - — = —^o, 119 
8 8m 

where we have inserted again ep{n) = (ir/m)n for the Fermi energy of free fermions at 
the carrier density n = n(T = 0). 

Note that this strong-coupling behavior of Tkt coincides roughly with the estimate of 
T c for three-dimensional superconductors in Refs. [[|, ||, and |J via the onset of Bose 
condensation of tightly bound, almost free composite bosons. In the first two of these 
references which include only quadratic fluctuations around the mean field (corresponding 
to ladder diagrams), T c was shown to tend to a constant value which does not depend 
on the internal structure of composite boson and is simply equal to the condensation 
temperature of a gas of free bosons of mass 2m and density n/2, implying that the 
interactions between the composite bosons is irrelevant in this approximation. 

We find the same situation in two dimensions, where Tkt tends to a constant depend- 
ing only on the mass 2m and the density n/2 of the pairs. No dependence on the coupling 
strength is left. The difference with respect to the three-dimensional case is that here the 
transition temperature T c = Tkt is linear in the carrier density n, while growing like n 2//3 
in three dimensions. Our result ( |119|) agrees with Ref. p6[ and . 



If interactions between condensed and noncondensed composite are taken into account, 
as done in Ref. || in three dimensions, then T c turns out to grow slowly with the coupling 
strength in the strong coupling regime. 



Equation ( |119| ) determines the critical temperature in the strong-coupling limit com- 



pletely in terms of the carrier density n . There exists a corresponding equation for the 
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temperature T* in the strong-coupling limit e 3> ep\ 

e 1 



(120) 



2 logeo/e^ 

For experimental purposes, the dependence of the ratio 2A(0) /T^t on the coupling 
strength is of interest. It is plotted in Fig. |T0|. Analytically, we have in the weak-coupling 
limit 

^}+(9(.r f7 3 ). (121) 



2A(0) 2n e 7 4 
-7jM = — i+ — + 



1 /4e 7 \ 2 



and for strong-couplings: 



2A(0) 
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-64x . 



(122) 



The two curves can easily be interpolated graphically to all coupling strengths, as seen in 
the figure. 
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Figure 10: Weak-coupling and strong-coupling estimates for the ratio 2A(0) /Tkt (solid 
curves). The dashed line is a graphical interpolation. 
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5 Phase Fluctuations in Three Dimensions and 
Superconductuctive Transition in 3D XY-Model 



In this chapter we discuss effects of fluctuations in three-dimensional systems. In three 
dimensions we have the expression for stiffness 

W , A) = ±n(,,T, A) - ^-jy-J&±^— ( 123 ) 

With it we can immediately set up 3D XY model, whose Hamiltonian reads: 

H = ^J d 3 x[W(x)] 2 . (124) 

Temperature of the phase transition of this model can be estimated using mean-field 
methods for the lattice 3D XY-Model §: 



T™ F ~ 3J 3D a, (125) 

a = 1/ n^ 3 is the lattice spacing of the theory || where is number of pairs. 
In the weak-coupling limit, the stiffness can be expanded near T* as follows: 



7 r> 3 A 2 

</ = ^tC(3)^ At, (126) 



This is similar to the coefficient of the gradient term in the Ginzburg-Landau expansion 
except that there it is obtained from a small-A expansion, whereas here the background 
gap has a nonzero modulus. 

Obviously we have no separation of the two temperatures in the BCS limit. For 
moderately strong coupling, the two temperatures are related by 

f c = f*- af* 5/2 , (127) 

where a = (2vr 2 ) 2 / 3 /3. 

In the strong-coupling limit of the theory where we have tightly bound composite 
bosons, the phase stiffness tends asymptotically to: 

V/i 2 + A 2 " 



= n_ _ 3^271-77^ 3/2 ex 
4m 167T 2 
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which may be expanded 

n 3 a /Onrrn _ ,_ V it.' 

(129) 



J= ^_3v^ r3/2exp 



4m 167T 2 
It obviously tends in this limit quickly to 



n 

J=—. 130 

4m v ' 

An estimate for the critical temperature, obtained via the mean-field treatment of the 3D 
XY-model on the lattice reads in this limit: 



2m 



^J-^T^expf-^ 
2 J n 1 '***/* c P l 



(131) 



This quickly tends from below to the value: 



*?™ = ^ = e W H0 - 19 *'- (132) 

This result is very close to the temperature of the condensation of bosons of mass 2m and 
density n/2, which, as it was discussed in the introduction was obtained including the 
effect of Gaussian fluctuations into the mean-field equation for the particle number [0, 
yielding Q 

T Bosons = [ n /2((3/2)] 2 / 3 7r/m = 0.218e F . (133) 

The XY-model nature of the phase transition at T c has been demonstrated in recent 
experiments on YBa2Cu307_«5 near the region of optimal doping. The phase transi- 



tion of this model was discussed in great detail on a lattice in the textbook || , since it 
describes the critical properties of the superfluid transition of Helium. 

6 Conclusion 

We have discussed the two different transitions taking place in superconductors at strong 
couplings, the formation of pairs and onset of phase coherence. For theory it was a 



6 When critical temperature is studied via retaining gaussian corrections to the number equation the 
crossover of the critical temperature has an artificial maximum in the region of intermediate couplings 
[0, El, so in this case it's limiting value is approached from above in the strong-coupling limit, this is not 
the case in our approach. 



33 



fortunate fact of history that the early-discovered metallic superconductors had such a 
weak coupling that there was only one transition which, moreover, can be understood by 
mean-field methods. In high-T c superconductors, the existence of a pseudogap brings in 
complications which we have tried to illuminate in the framework of the simple fermion 
model with 5-function interaction. The interpretation of the experimental data is still 
complicated due to the complex chemical structure of these compounds. Little is known 
up to now on the real forces causing the pairing. 

We have studied the crossover from BCS-type to Bose-type superconductivity and 
the behavior of the pseudogap state. Our crossover parameter is xq, a quantity closely 
related to the chemical potential. For this purpose we have used the gradient expansion 
of the effective energy functional to set up an equivalent XY-model which allows us to 
investigate the onset of long-range order in the phase fluctuations. In two dimensions, we 
have given a simple analytic expression which shows how the resulting Kosterlitz-Thouless 
temperature Tkt at which quasi-long-range order sets in moves towards the pair-binding 
temperature T*, and merges with it in the weak-coupling limit. We have found similar 
results in three dimensions, setting up a three dimensional XF-model for the description 
of the onset of phase coherence in the superconductive transition, and how this transition 
evolves to the ordinary BCS transition in the weak-coupling regime. 

We have also studied the weak- to strong-coupling crossover of thermodynamic func- 
tions of the superconductive state near zero temperature as well in the pseudogap phase 
near the critical temperature. 

Certainly, our mean-field estimates for T* are quite crude, and we expect significant 
modifications due to fluctuations, in particular of the character of the transition which 
experimentally does not seem to be of second order, and may not be a phase transition 
after all. All our formulas for thermodynamic quantities will have to be smeared out in 
temperature near T*, before any possible comparison with experiments. 

Let us finally remark that the separation of T* and T c has an analogy in the ferro- 
electrics and magnets which also contain two separate characteristic temperatures, for 
example in the later case — the Stoner- and the Curie-temperature. 
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A Action functional of Collective Pair Field 

In this appendix we briefly outline derivation of the effective action (P5j). As shown in 



Ref. |p3[ , a pair field A is introduced to eliminate the quartic interaction term in the 
functional integral involving the action of the Hamiltonian (JJ): 



A= dt 



[d D x^i(x)thd t Mx)-H(t) 

ST J 



(134) 



After that, the fermions can be integrated out. At a constant pair field, we find the 
potential part (PP| ) of the collective-field action 

f IAI 2 + °° r d D k 

n pot (»,T,A,A*) = V IL-L-T J2 J T^ptrllnG- 1 ^,^^^ 3 ] 



9 n ^oo J ( 27r ) 

dPh 
W) 1 



+ T E J 7^y5^ogG^^ n , k)e**"»] , 5 - +0, (135) 



where A = Ae td and 



G~\iu n , k) = iuj - r 3 £(k) + r + A + r_A* = ( *"» ^ (k) ) . (136) 

Using the identity trlogA = logdet A, equation ( |135| ) takes the form 

O ( T A AM T//I A ! 2 T V 5 f 1 ^ + ^(k) + |A| 2 

n pot (fi,T,A,A = v|— - - r^y—tog—--— 

" /(Cp K(k) + £(k)] } ! (137) 
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After performing the sum over the Matsubara frequencies in (|137[) , we obtain the well- 
known mean-field expression for f2 pot 



n pab ( f i,T,A,A*) = V 



+ 



v 9 
d D k 



d D k 



£ 2 (k) + |A| 2 
2T log 2 cosh v ^ - f (k) 



2T 



2Tlog2cosh^-e(k) 



(138) 



(2n) D 

In order to derive the kinetic part f2kin of the mean-field energy, we must calculate the 



first two terms of the series (|98[), the first being |24] 



7-D, 



T 



(27T) 



with 



G(iui n , k) 



£ d D ktT[g(icu n ,k)r 3 

=-oo 

iu; n I + r 3 £(k) - nA 



ia r g + (w) 2 



^ 2 +e 2 (k) + A 2 • 

After summing over the Matsubara frequencies and integration over k, we obtain 



fig) = T jf drj d D xn(fi, T, A) 



zd T fl | (W) 2 



5m 



(139) 



(140) 



(141) 



with n(/x, T, A) given in by (§). After an ansatz £ = r 3 Oi + IO2, where Oi and O2 are the 
two gradient terms in Eq. (|101|) , we derive for the second term the two contributions 
from 0\ and 2 : 



nS(Ox) = / d D x K(fi,T, A) 



id T # (V0) ; 



5m 



(142) 



where K(fi,T, A) was given in ( Jl 04 ) . This is the second term in ( |102| ). The second term 
in ( |103| ) is obtained from the second contribution to fij^: 

- f dr [d D x—^ f ^ k " 
Jo J 32vr 2 m 2 



d u k- 



cosh 2 [J£ 2 (k) + A 2 /2T] 



(V9f 



(143) 



Combining (Q, (|T4|) and (HJ) we obtain ([1021) . 
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